Doping-induced vertical line nodes in the superconducting gap of the iron arsenide 
Bai ^.K^Fe2As2 from directional thermal conductivity 
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The thermal conductivity k of the iron-arsenide superconductor Bai_a:KiFe2As2 was measured 
down to 50 mK in a magnetic field up to 15 T, for a heat current paraUel and perpendicular to 
the tetragonal c axis. In the range from optimal doping {x ~ 0.4) down to a; = 0.16, there is no 
residual linear term in ^(r) as T — >■ 0, showing that there are no nodes in the superconducting gap 
anywhere on the Fermi surface. Upon crossing below x — 0.16, a large residual linear term suddenly 
appears, signaling the onset of nodes in the superconducting gap, most likely vertical line nodes 
running along the c axis. We discuss two scenarios: (1) accidental nodes in an s-wave gap, resulting 
from a strong modulation of the gap around the Fermi surface, in which minima deepen rapidly 
with underdoping; (2) a phase transition from a nodeless s-wave state to a d-wave state, in which 
nodes are imposed by symmetry. 

PACS numbers: 74.25.F-, 74.20.Rp 



To elucidate the pairing mechanism in iron-based su- 
perconductors, a key question is whether the underlying 
interaction is anisotropic. This will typically show up 
as anisotropy in the gap, as in cuprate superconductors 
whose d-wave gap goes to zero along lines of nodes im- 
posed by symmetry. The superconducting gap of iron 
pnictides has been investigated with several probes, but 
a clear overall picture has yet to emerge [T]. In some 
materials, like LaFePO 0, KFe2As2 [S] and P-doped 
BaFe2As2 [4], nodes have been detected. In other ma- 
terials, like LiFeAs E] and Co-doped BaFe2As2 (Co- 
Bal22) [7], an isotropic gap without nodes is seen. How- 
ever, the gap in Co-Bal22 is only isotropic at optimal 
doping, where Tc is maximal, and it acquires nodes when 
doping is either reduced or increased [8, 9 . While several 
calculations favour a nodeless sign-changing s± state, the 
gap can have strong modulations and even nodes [TnHT5] , 
and a d-wave state has comparable energy [TTl [T^ ■ 

In this Letter, we investigate the gap structure of K- 
doped BaFe2As2 using low-temperature thermal conduc- 
tivity, a bulk directional probe of the superconducting 
gap |14j . Our main finding is the sudden appearance of 
a finite residual linear term, kq/T, upon crossing below 
the critical concentration x — 0.16, the signature that 
nodes have appeared in the superconducting gap. We 
discuss whether these are accidental, the result of deep 
minima in a strongly-modulated s-wave gap, or imposed 
by symmetry in a doping-tuned transition to a d-wave 
state. In either case, our findings reveal how easily the 
full isotropic gap observed at maximal Tc can be desta- 
bilized by tuning away from optimal conditions, pointing 



to anisotropic and competing interactions. 

Experimental. - Single crystals of Bai_2,Kj;Fe2As2 (K- 
Bal22) were grown from self flux [TB]. The doping level 
X was determined from the superconducting transition 
temperature Tc, defined by the point of zero resistance, 
using an empirical curve of Tc vs x [13 [T5] (see Fig. 1). 
Tc values for a-axis samples are: 9.4, 10.0, 14.5, 16.4, 
25.5, 29.5, 36.7 K; Tc values for c-axis samples are: 10.7, 
11.3, 17.1, 18.8, 28.6, 28.9, 38.1 K. Our samples cover 
a wide range, both inside and outside the region of co- 
existence with antiferromagnetic order (below x ~ 0.24 
|16)). Thermal conductivity was measured in a standard 
one-heater two-thermometer technique for two directions 
of heat flow: parallel (J||c; Kc) and perpendicular (J||a; 
Ka) to the [001] tetragonal c axis. Sample geometry, con- 
tact technique and measurement protocol are described 
elsewhere [S]. The magnetic field H was applied along 
the c axis. 

Results: H — 0.- The thermal conductivity of K- 
Bal22 is displayed in Fig. 2. The fits show that the 
data below 0.3 K are well described by the function 
k/T = a + bT°'. The first term, a = kq/T, is the resid- 
ual linear term, entirely due to electronic excitations |14j . 
The second term is due to phonons, which at low tem- 
perature are scattered by the sample boundaries, with 
1 < Q < 2 [mill]. For X > 0.16 {Tc > 13 K), kq/T = 
within error bars, for both current directions. This is 
consistent with our previous data for Ka at Tc=26 K 
(Ref. . By contrast, for both directions a large resid- 
ual linear term is found when x < 0.16 (Tc ~ 10 K). 

It is convenient to compare kq/T to the normal-state 
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FIG. 1: Top: Doping phase diagram of K-Bal22, showing 
the antiferromagnetic (AF) phase below Tn (blue squares, 
from ref. ITS)) and the superconducting (SC) phase below 
(open (closed) red circles, for o-axis (c-axis) samples), which 
coexist in the shaded region. Bottom: Residual linear term 
in the thermal conductivity k as T — >■ 0, no/T, plotted as a 
fraction of the normal-state conductivity, kn /T, for both Ka 
(open symbols) and Kc (closed symbols), for magnetic fields 
H=0 (red) and = 0.15 Hc2 (black). Error bars reflect the 
combined uncertainties in the extrapolation of k/T and p to 
T = 0. The vertical dashed line at a; = 0.16 separates a region 
in which the superconducting gap has nodes [x < 0.16) from 
one in which it doesn't (x > 0.16). The solid black line is a 
guide to the eye. There is no superconductivity in the shaded 
region. 



conductivity kj^/T, estimated using the Wiedemann- 
Franz law, namely kn/T — Lq/pq where Lq = 
(7r^/3)(fcB/e)^, applied to the extrapolated residual re- 
sistivity pq, as discussed in Ref. HI In Fig. 1, the ra- 
tio (kq/T)I{k^/T) is plotted vs x. It jumps suddenly 
from to ^ 0.5 upon crossing below x ~ 0.16. This 
signals the onset of nodes in the superconducting gap of 
K-Bal22. The fact that {kq/T) / [k^ /T) is of the same 
magnitude for both current directions indicates that the 
nodes preserve the a — c anisotropy of normal-state trans- 
port, consistent with line nodes that run vertically along 
the quasi-2D Fermi surface of K-Bal22. 

Results: H > 0.- In Fig. [3j we plot the residual linear 
term kq/T as a function of magnetic field H. As in Fig. 1, 
Ko/T is normalized to kn/T, and H is normalized to Hc2, 
the upper critical field at T — )■ 0, given approximately by 
= 100 T {Tc/AO K) [ini m]. At high x, the field 
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FIG. 2: Thermal conductivity k{T) of K-Bal22 at three rep- 
resentative K concentrations, indicated by their Tc values, 
plotted as k/T vs T for a current parallel (J||c, right panels) 
and perpendicular (J||a, left panels) to the tretragonal c axis. 
Data are shown for three values of the magnetic field, as indi- 
cated (data for other fields are not shown for clarity). Lines 
are a fit to ft/T — a + bT" , used to extract the residual linear 
term a = kq/T, plotted vs x in Fig. 1 and vs H in Fig. 3. 



dependence is very weak, either flat at Tc = 38 K or with 
a slight upward curvature at Tc — 30 K. This is typical 
of s-wave superconductors, in which a full gap is present 
on the entire Fermi surface. Quasiparticle transport is 
caused by tunneling between states localized in adjacent 
vortex cores and thus grows exponentially as the intervor- 
tex separation shrinks with increasing field [22 . By con- 
trast, when X < 0.16, kq/T increases rapidly at low H, 
with the downward curvature characteristic of a d-wave 
superconductor such as the overdoped cuprate Tl-2201 
P5] (Fig. 3). In this case, delocalized nodal quasiparti- 
cles outside the vortex cores are excited by the Doppler 
shift in their energies caused by the field [24j. 

In the range above x = 0.16, throughout which the gap 
has no nodes, the field dependence varies significantly, 
from flat at a; ~ 0.4 to strong at a; ~ 0.18. In samples 
with Tc ~ 15 K, the H dependence is very similar to that 
of a multi-band s-wave superconductor like NbSe2 [22] 
(Fig. 3) , in which the small value of the superconducting 
gap on one part of the Fermi surface leads to enhanced 
excitation. In Fig. 1, we plot the x dependence of kq/T 
measured at H = 0.15 Hc2- We see that the sudden 
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FIG. 3: Residual linear term kq/T in K-Bal22 normalized 
by the normal-state conductivity, kn/T, as a function of mag- 
netic field H normalized by the upper critical field Hc2- The 
data for in-plane (J||a, open symbols) and inter-plane (J||c, 
closed symbols) transport are shown for five K concentrations, 
indicated by their Tc values. For comparison, we reproduce 
corresponding data for the isotropic s-wave superconductor 
Nb (see ref. lT4|l . the multi-band s-wave superconductor NbSe2 
[22| . and the d-wave superconductor Tl-2201 [23] . 



jump in kq/T at = 0, when only zero-energy nodal 
quasiparticles contribute, turns into a more gradual rise 
at finite iJ, when finite-energy quasiparticles are excited 
across a small gap. Note that there could still be a small 
jump across x = 0.16 in finite H . 

Coexisting AF order.- Neutron studies show that an- 
tiferromagnetic (AF) order is present down to at least 
T ~ 1 K from a; = up to at least x ~ 0.24 (T^ = 26 K), 
and it occupies at least 95 % of the sample volume (TB]. 
The fact that kq/T = (at = 0) for > 13 K im- 
plies that putative normal-state regions (with kq/T = 
k-m/T > 0) must occupy a volume fraction below the 
percolation threshold. These two constraints rule out a 
scenario of inhomogenous phase separation in the coexis- 
tence region, whereby the sample would consist of mutu- 
ally exclusive regions of pure AF order and pure super- 
conductivity. Therefore, there is a substantial region of 
the phase diagram where nodeless superconductivity co- 
exists homogeneously with AF order. It has been argued 
that this is consistent with a s± state but inconsistent 
with a standard s++ state . 

Scenario I: extended s-wave.- Given that the field- 
induced enhancement of kq/T increases as x approaches 
the critical value {x = 0.16) below which nodes appear 
in the gap (Fig. 1), it is natural to view the in-field be- 
haviour as a precursor to the nodes, namely as minima 
in the gap that become gradually deeper until they touch 
the Fermi surface. This is the so-called 'extended s-wave' 
scenario, illustrated in Fig. 4. 
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FIG. 4: Two scenarios for the doping evolution of the super- 
conducting gap structure in K-Bal22 (solid red line), sketched 
for a single 2D Fermi surface (dashed circle) . In the extended 
s-wave scenario (top): a) the gap is isotropic at high x; b) 
it develops minima in certain directions at intermediate x; c) 
those minima deepen to become negative at low x, producing 
nodes where A = at 8 points on the Fermi surface. In the 
d-wave scenario (bottom), a phase transition occurs between 
a superconducting state with s-wave symmetry at high x (d) 
and one with d-wave symmetry, and associated nodes, at low 
X (e). 



The isotropic gap at high x develops minima in certain 
directions as x decreases, until this modulation as a func- 
tion of azimuthal angle becomes strong enough that the 
minima become negative, thereby producing nodes below 
some critical x value. These nodes are called 'accidental', 
because they are not imposed by symmetry. There is no 
transition across the critical doping where nodes appear - 
the entire evolution takes place within one phase, with no 
change of symmetry. The fact that the in-field behavior 
at a; > 0.16 is just as isotropic as the zero-field behaviour 
at X < 0.16 reinforces the connection between minima 
and nodes. It is possible that the reconstruction of the 
Fermi surface caused by the onset of AF order (below 
X ~ 0.24) either triggers or favours the gap modulation 
we observe, first as minima then as nodes. 

In a number of calculations applied to pnictides, a 
superconducting state with s± symmetry is found to 
be the most stable, but often the associated gap func- 
tion has strong modulations, possibly leading to acci- 
dental nodes [TTtiT5] . The gap modulation comes from 
a strongly anisotropic pairing interaction, which is also 
band-dependent. In 2D models, it is typically the gap on 
the electron Fermi surface centered at the M point of the 
Brillouin zone which shows a strong angular dependence 
within the basal plane [11] [T^], leading to what would 
be vertical line nodes in 3D. In 3D calculations, nodes 
have been found on the hole Fermi surface at the F point 
[27] . This raises the question of which sheet of the Fermi 
surface of K-Bal22 has gap nodes. (Note that the nodes 
occur inside the phase of coexisting AF order, on a Fermi 
surface that is reconstructed relative to the paramagnetic 
Fermi surface assumed in the calculations.) 
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Because nodal conduction is ~ 50 % of the normal- 
state conductivity (Fig. 3), which is dominated by hole- 
like carriers (Hall [28] and Seebeck [29] coefficients are 
positive), the dominant hole-like sheets must have nodes 
(and possibly other sheets as well). 

Comparison with Co-Bal22.- There are clear similar- 
ities in the behavior of kq/T in Co-Bal22 and K-Bal22. 
At optimal doping, the gap is nodeless in both cases, 
with kq/T = Q in both directions. The field dependence 
is more rapid in Co-Bal22, indicative of gap minima al- 
ready present at maximal Tc, while they are absent in 
optimal K-Bal22. This points to a correlation between 
gap modulation and Tc weakness, as T™^'^ = 25 K in 
Co-Bal22 and T™'*'^ = 38 K in K-Bal22. Away from 
optimal X, the gap minima deepen for both materials. In 
Co-Bal22, this is seen not only on the underdoped side 
(with the coexisting AF order) , but also on the overdoped 
side, where there is no AF order. In both cases, the 
gap minima lead to isotropic in-ficld transport (relative 
to kn)j and so are consistent with vertical line minima 
along the c axis. 

The difference lies in the nodal character. While K- 
Bal22 has vertical line nodes, Co-Bal22 only has nodes 
on regions of the Fermi surface that contribute strongly 
to c-axis conduction. Recent calculations can account 
for this in terms of a Fermi surface with strong c-axis 
dispersion [3D]. In a scenario where nodes are the result 
of minima that gradually deepen to become negative, this 
deepening would be uniform in K-Bal22, independent of 
kz location on the Fermi surface, but fc^-dependent in Co- 
Bal22. It remains to be seen whether the nodes detected 
in KFe2As2 via a large in-plane kq/T [3] can be viewed 
as the end point of a gradual deepening of line minima 
in the overdoped regime of K-Bal22. 

Scenario II: transition from s-wave to d-wave.- In a 
second scenario, the nodes appear below x — 0.16 be- 
cause of a phase transition from an s-wave state (e.g. 
s±) to a d-wave state. The nodes would no longer be ac- 
cidental, but imposed by the d-wave symmetry, and they 
would necessarily appear on the Fermi-surface sheets cen- 
tered on r, presumably hole-like. This would be compat- 
ible with the large value of kq/T and the dominance of 
hole-like carriers in the normal-state conductivity. The 
strong field dependence of kq/T for x above the would- 
be transition a,t x = 0.16 requires that the s-wave gap 
develops deep minima before the transition to d-wave. 

In summary, the thermal conductivity of K-Bal22 in 
the r = limit reveals that nodes appear in the su- 
perconducting gap when the K concentration is reduced 
below a critical value. Because the nodal quasiparticles 
are as effective in conducting heat parallel and perpendic- 



ular to the c axis, the nodes are likely to be vertical lines 
on the quasi-2D Fermi surface, running parallel to the c 
axis. The large value of the nodal conductivity implies 
that nodes are present on those hole-like sheets of the 
Fermi surface that dominate the normal-state transport. 
The fact that the isotropic gap observed in K-Bal22 when 
Tc is maximal becomes anisotropic as superconductivity 
weakens suggests that pairing in iron-arsenide supercon- 
ductors involves competing interactions. 
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